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Summary 

To easily study combustor design parameters using computational fluid dynamics codes (CFD), a 
Gridgen Glyph-based macro (based on the Tel scripting language) dubbed BladeMaker has been 
developed for the meshing of an idealized, single radial swirler, lean-direct-injection (LDI) combustor. 
BladeMaker is capable of taking in a number of parameters, such as blade width, blade tilt with respect to 
the perpendicular, swirler cup radius, and grid densities, and producing a three-dimensional meshed radial 
swirler with a can-annular (canned) combustor. This complex script produces a data format suitable for 
but not specific to the National Combustion Code (NCC), a state-of-the-art CFD code developed for 
reacting flow processes. 


Introduction 

Lean direct injection (LDI) is a NASA-developed, low-emissions gas turbine combustion concept. 
Although this concept produces reduced emissions in flame tube tests, the mixing and combustion 
processes must be understood to optimize the concept (refs. 1 to 3). Computational fluid dynamics (CFD) 
programs allow LDI design studies on computers. However, in past studies, the geometries were fixed 
because of the difficulties in creating a computational mesh (refs. 4 and 5). To better optimize the LDI 
concept using current CFD codes, parametric geometry manipulation is a necessary requirement. 

Although parametric geometries may be easily created by using computer-aided design (CAD) tools 
like Pro Engineer (Parametric Technology Corp.), creating a suitable efficient mesh for CFD design 
studies is difficult. Coupling parametric solid models to a mesh generator and creating a good quality 
mesh are still a research problem. To solve this problem, a Gridgen (Pointwise, Inc., ref. 6) Glyph macro 
dubbed BladeMaker has been developed for the meshing of a simplified, idealized, LDI radial swirler. 
BladeMaker is capable of taking in a number of parameters, such as blade width, blade tilt with respect to 
the perpendicular, swirler cup radius, and grid densities, and producing a meshed radial swirler. This LDI 
mesh can then be used as input to a CFD code, preferably the National Combustion Code (NCC). 

Briefly, The NCC is a state-of-the art CFD program specifically designed for combustion processes. 
A short summary of the features of NCC follow: the use of unstructured grids (ref. 7); massively parallel 
computing — with almost perfectly linear scalability (refs. 8 and 9); a dynamic wall function with the 
effect of adverse pressure gradient (ref. 10); low-Reynolds-number wall treatment (ref. 1 1); a cubic, 
nonlinear k-epsilon turbulence model (refs. 12 and 13); and stiff laminar chemistry integration. Recently, 
viscous low-speed preconditioning (refs. 14 and 15) has been added to improve the low-speed 
convergence of the NCC in viscous regions. The combination of these features is usually not available in 
other CFD codes and gives the NCC an advantage when one is computing turbulent, reacting flows. 
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NOTES ON FORMATTING AND NAMING: This document uses different fonts to indicate 
specific code items. Sample code will always be shown in 12-point Courier New font face and 
will be boxed. Commands to be submitted to a shell and file names will be in unboxed 1 2 -point 
Courier New . Inline Glyph code or commands will also be in this font, but not boxed. Glyph 
variables and procedure names will be shown in 1 0-point Arial when not in sample code. 

Using BladeMaker 

Creating a Project Directory 

It is suggested that you create a directory for your case inside the BladeMaker directory. This 
suggestion is made because modification of the script and input files will be required. If you begin to have 
problems with a file, you will always have the original to go back to. Once you have created a project 
directory, you will want to copy the files param . in and all the files with the extension sh into it. 

These files can be found in the BladeMaker root directory. 

Parameter File Construction 

All specifications for your geometry and mesh are made in the parameters file, which is called param . in 
by default. The parameter file that comes with BladeMaker is shown in figure 1 for quick reference. 

ASW_TYPE 

NCC 

ASW_DIM 

3 

swirler_rad 

2.0 

bladevol_rad 

6.0 

manifold_rad 

9.0 

num_b lades 
8 

blade_ang 

15.0 

blade_rad 

0.1 

height_s 

1 . 0 

basic_con_dim 

5 

blade_con_dim 

10 

wall_ds 

O.OOl 

height_cc 

15.0 

fname 

/ your /working/ di rectory /pa tr an . out 
*EOF* 

Figure 1. — Default parameter file, param. in. 
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TABLE 1.— DESCRIPTIONS AND LIMITS IN PARAMETER FILE 


Variable name 

Description 

Limit 

Data type 

Lower | Upper 

ASW TYPE 

Analysis software package 

Any valid Gridgen analysis software type. See Gridgen 
documentation for more infonnation. 

ASW DIM 

Number of dimensions in 
geometry 

3 

3 

Positive integer 

swirler rad 

Radius of swirler cup 

1.0 

Less than 
bladevol rad 

Positive real 

bladevol_rad 

Distance from center of 
swirler cup to edge of 
blade volume 

Greater than 
swirler rad 

Less than 
manifold rad 

Positive real 

manifold rad 

Radius of entire radial 
swirler 

Greater than bladevol rad 

Infinite 

Positive real 

num blades 

Number of blades in swirler 

These two variables are inversely proportional. The minimum number 
of blades possible for any blade rad is 3. The macro will fail to run on 
two blades for unknown reasons. 

blade_rad 

One-half of blade width 

blade_ang 

Tilt of a single blade with 
respect to the perpendicular 
(Counterclockwise is 
positive and clockwise is 
negative.) 

-50.0 

50.0 

Real 

height_s 

Height of region of swirler 
containing the blades 

0.0 

Infinite 

Positive real 

height_cc 

Height of combustion 
chamber 

0.0 

Infinite 

Positive real 

basic con dim 

Number of points placed 
initially on a connector 

2 

Infinite 

Positive integer 

blade_con_dim 

Number of points placed on 
a connector near blades 

Greater than 
basic con dim 

Infinite 

Positive integer 

wall ds 

Average spacing for grid 
points near walls 

Greater than 0.0 

0.01 

Positive real 

fname 

File name for analysis 
software file (i.e., patran.out) 

Any writable absolute file name 

*EOF* 

This is not really a variable, but it signals the end of the file. | 


The top section of the file par am . in serves as an identification of its contents and is ignored by 
BladeMaker. Each line after that is either a variable name or a value. If a variable name is encountered, 
the parser sets it to the value on the next line. So, for the file in figure 1, ASW_TYPE is set to NCC, 

ASW DIM is set to 3, and so on. A description of each of these variables and the limits of BladeMaker are 
found in table 1. Because of “stretching” of the structured mesh, the variable blade_ang has a limitation 
of approximately -50 to 50° to the perpendicular to keep mesh quality high. 

Glyph Script 

Most users will only have to touch one line in the actual Glyph script, blademaker . gif. This is 
a plain text file and can be edited in any text editor. Just one line has to be modified in this script, as 
highlighted in figure 2. PARAM FILE must be set to the absolute file path of your parameter file instead 
of to the default / your /working/ directory/ param . in. 
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######################################## 

# PARAMETERS # 

######################################## 

# Set Pi - Tel uses radians in trig functions, 
set PI [expr {4.0 * atan (1.0)}] 


set PARAM_FILE " /your/working/directory/param . in 
set EOF " *EOF* " 


Figure 2. — Changes in the glyph script blademaker . gl f . 


TABLE 2.— SHELL SCRIPT DESCRIPTIONS 


Shell script 

Description 

Output 

unix-unopt . sh 

Runs BladeMaker 

An analysis software boundary conditions file, such 
as patran . out 


# ! /bin/bash 

# BladeMaker 

# Gridgen Only Run Shell Script 

# UNIX Version 

# Created 19 July 2005 

source /usr/local/etc/setup . sh 
setup gridgenvl57 
gridgen blademaker.gif 


Figure 3. — The unix-unopt . sh bash shell script. 


Running BladeMaker 

BladeMaker has been run as a batch script since version 1.0.2. We assume that Gridgen Version 
15.07 or later is available on your computer. Simply instruct your shell to interpret the shell script 
appropriate for your case. Table 2 gives a description of the shell scripts. 

The shell script unix-unopt . sh (fig. 3) will run BladeMaker alone. It does not li nk 
automatically to CFD codes or optimizers. This script assumes that you have access to the bash shell and 
have available Gridgen Version 15.07 or later. If your system is not configured as such, you will need to 
use this section as a basis for creating your own shell script. 

The first line sources in the setup . sh shell script. On the second line, setup . sh is used to 
make Gridgen Version 15.07 available. The third and final line of code runs Gridgen in batch mode on the 
macro in blademaker . gl f . A boundary condition file useable by an analysis software package will 
be created. This is usually patran . out, a PATRAN neutral file useable by the NCC. 
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Breakdown of BladeMaker Glyph Script 

Previous Knowledge 

The appendix assumes that the reader has basic knowledge of Glyph, or at least Tel. Some sections 
may be confusing without such knowledge. More information on Tel can be found at http://www.tcl.tk/ 
and more information on Glyph can be found in Gridgen Version 15: Glyph Reference Manual. 

General Overview 

The BladeMaker script itself is stored in the ASCII text file blademaker .gif. Because the 
contents of blademaker . gif are referenced throughout the report, blademaker .gif is listed in 
the appendix. It is divided into a number of sections, indicated by comment headers. These sections 
appear in the following order: parameters, procedure definitions, and the main Glyph script. Each is 
covered in detail in the next sections. 


Parameters Section 

The purpose of the parameters section is twofold. First, three constants referred to by the macro are 
defined. Since Tel does not provide a built-in re, the constant PI is set to 4.0 * arc tan 1.0, which is a good 
approximation of 7t. RUN MODE defines whether or not the Gridgen will use the script interactively, set 
to 1, or in a batch mode, set to 0. PARAM FILE once again defines the input parameters. 

Procedure Definitions 

In the procedure definitions section, all procedures used in the macro are defined. The section is 
subdivided into utility, connector, database, meshing, boundary condition, mesh refinement, and large- 
scale creation procedures. 

Utility procedures . — Utility procedures are called in the procedure calls section but do not actually 
create any grid structure. Currently, the only procedure in this section is ClearWorkspace. This procedure 
takes no arguments and prepares the workspace for use by the macro by clearing any grid structure and 
setting all options to their default values. It also sets the analysis software type and dimension as defined 
in the parameters file. 

Connector procedures . — Connector procedures bundle groups of primitive Glyph commands for 
certain functions. Without this section, these groups of commands would have to be repeatedly rewritten, 
which is tedious and leads to a large file size. These procedures are usually only called from within the 
large-scale creation procedures. 

The first procedure in this section, SplitConlntoNEqualParts, was adapted from the Glyph procedure 
ConSplitlntoN found on the Pointwise Glyph exchange at http://www.pointwise.com. The procedure 
SplitConlntoNEqualParts takes a positive integer as the first parameter and a list of Glyph connector IDs 
as the second parameter and splits each connector in the list into the specified number of equal parts. It 
also returns the connector IDs of the resultant connectors. 

The next connector procedure is CreateEndptCon. It takes a list of positive integers as its only 
parameter, which must have an even-numbered length. For each successive pair of positive integers in the 
argument, a connector is created between the endpoints of the corresponding elements in the list. This is 
returned by the Glyph command conGetAll in the namespace gg. For example, if we passed this 
procedure the list 0 3 10 63, two connectors would be created. The first would be created between the 
endpoints of the connectors with Glyph IDs [lindex [gg :: conGetAll ] 0] and [lindex 

[ gg : : conGe tAl 1 ] 3 ] , and the second connector would be formed between the endpoints of the 
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connectors with Glyph IDs [lindex [gg : : conGetAll ] 10] and [lindex 

[gg : : conGetAll ] 63]. 

The final connector procedure is Create2PtCon. This procedure takes as a parameter a list of points in 
three-dimensional space, which must have an even- valued length. A straight connector is drawn between 
the first point and the second point, then the third point and the fourth point, and so on. For example, if it 
were passed the list [[0 4 7] [8 19 23] [7 5 8] [0 0 0]], two connectors would be 

created. The first would start at (0,4,7) and end at (8,19,23). The second would begin at (7,5,8) and end at 
the origin. 

Database procedures . — These procedures are much like the connector procedures, except that they 
are used to create database entities instead of connectors. 

Currently there is only one database procedure called Create2PtDB, which takes as a parameter a list 
of indexes to the return value of dbGe t All. This parameter must have an even-numbered length, and 
these indexes must correspond to the Glyph IDs of points. The list is treated much like the list in 
Create2PtCon, except that straight database entities are created instead of straight connectors. For each 
pair in the parameter, a straight database entity is created beginning at the first element in the pair and 
ending at the second element. For example, if we passed this procedure the list [0 1 34 68], two 
straight database entities would be created. The first would be between the points with Glyph IDs 0 and 1, 
and the second between points with Glyph IDs 34 and 68. 

Meshing procedures . — This collection of procedures serves to bundle the primitive commands used 
to create structured mesh entities such as structured domains into single procedures that can be easily 
called from the large-scale mesh creation procedures. 

The first such procedure is called CreateStrucDom and is used to create structured domains. The 
parameter for this procedure is a list of indexes to the return value of conGe t A1 1 . The parameter must 
have a length that is evenly divisible by 4. Each set of four indexes is treated as a definition of the edges 
of a single structured domain. For example, if this procedure were passed the list [0 4 13 55 0 8 
16 3 ] as input, two structured domains would be created. For the first domain, the first edge would be 
the connector with the Glyph ID found in the 0th element of the return value of conGetAll. The 
Glyph ID for the second edge would be found in the 4th element, the third in the 1 3th element, and the 
final edge in the 55th element. The second domain would be constructed in much the same fashion. This 
procedure eliminates the tedious copying and pasting of the same series of Glyph commands. 

The next and final procedure is CreateStrucBIk. This is identical to CreateStrucDom, except that it 
creates structured blocks and thus must take a list that is evenly divisible by 6 as input. The list elements 
are treated as indexes to the return value of the command domGetAll. 

Boundary condition procedures . — These procedures are for tagging boundary conditions as well as 
for generating the analysis software file, like patran . out. 

The first such procedure, called TagBoundaryConditions, tags certain domains as being boundary 
conditions. These boundary conditions are very specific to this particular test case. The bottom of the 
swirler (i.e., the lower domain in the blade region) is marked as having an “inlet-uvw” boundary. The top 
most domain on the combustion chamber that runs parallel to this is tagged as the “exit” BC. From a top- 
down view, the leftmost domain in the swirler cup is tagged as a “periodic-a” BC and the rightmost is 
tagged as a “periodic-b” BC. The leftmost blade volume domain is tagged as a Gridgen generic boundary 
condition “BndCond 20,” the leftmost combustion chamber domain as “BndCond 21,” the rightmost 
blade volume domain as “BndCond 22,” and finally the rightmost combustion chamber domain as 
“BndCond 23.” 

The final boundary condition procedure, called ExportASW, takes a file name as a parameter and 
creates the analysis software boundary condition file. If our analysis software is set to NCC, this would be 
the patran . out file. The file name tells where to create the file. 
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Mesh refinement procedures . — The procedures are for refining the mesh that has been previously 
generated. The first such procedure, called ClusterPts, is called from within other mesh refinement 
procedures. It takes three parameters, the average As value near the walls (ds), a list of indexes to the 
return value of conGe t A1 1 (index), and an end of the connector to cluster points near (loc). For each 
element of the return value of conGe t A1 1 that is listed in index, the procedure cluster points to a value 
of the average As value ds near the end specified in loc — which is either “start” or “end.” 

RefineSwirlerGrid takes two parameters, the average As value near the walls (wall_ds) and a connector 
dimension (blade_con_dim). This procedure increases the grid density near the blades in the lower 
portion of the swirler by changing the connector dimension from the default to blade_con_dim and cluster 
points near the blade walls by setting the average As value of the subconnectors pointing to the blade 
walls to wall_ds. 

RefineCombustionChamberGrid takes a single parameter, an average As value (wall_ds) and increases 
the density of the combustion chamber grid by a factor of one-quarter times the ratio of the height of the 
combustion chamber to the height of the swirler cup. This is done since the combustion chamber is 
usually much taller than the swirler cup, and the grid would otherwise become less dense in the 
combustion chamber. Next, points are clustered near the combustion chamber outside wall and the wall of 
the combustion chamber facing the swirler cup by setting the average As value of the appropriate 
subconnectors to wall_ds. 

Points must also be clustered along the wall in the buffer region between the swirler cup and the 
combustion chamber. There is one such wall facing the combustion chamber and another facing the 
swirler cup. The reader may inquire as to why there is no such procedure in blademaker . gif. The 
reason is that the procedure CreateBuffer assumes that both CreateSwirler and RefineSwirlerGrid have 
already been called. Thus, CreateBuffer creates a buffer volume that already has a refined mesh. 

Large-scale creation procedures . — These procedures call the database, connector, mesh creation, and 
mesh refinement, in addition to other smaller procedures, in such an order as to create large sections of 
the final grid and geometry. Although many of the previous procedures could easily be transplanted to 
another Glyph mesh generation script for a completely different case, these are specific to this geometry 
and mesh. 

The first of these procedures, MakeWhole, is an optional procedure in place for demonstration 
purposes only. As distributed, blademaker . gif never actually calls this procedure. However, the 
user could easily provide such a call by adding the line MakeWhole to the procedure calls section. This 
procedure simply creates the whole radial swirler from the slice via copying by rotation. 

The next procedure, CreateSwirler, is responsible for creating the lower portion of the radial swirler 
slice that contains the swirler cup and blades. It takes as parameters the radius of the swirler cup, the 
distance from the center of the swirler cup to a point on the edge of the blade volume, the radius of the 
whole radial swirler, the number of blades in the swirler, the angle a single blade makes with respect to 
the perpendicular, half the width of a single blade, the height of the swirler cup, and the number of points 
to initially place on all connectors in this region. The blade angle is considered positive in the 
counterclockwise direction and negative in the clockwise direction. 

The third procedure, CreateBuffer, takes just one parameter, the height of the swirler cup. It creates a 
slice of the buffer region between the combustion chamber and the swirler cup by copying the swirler cup 
via a translation of height equal to the parameter. 

The final large-scale creation procedure, CreateCombustionChamber, is responsible for the creation 
of the combustion chamber slice. It takes three parameters: the height of the swirler cup, the desired 
height of the combustion chamber, and the number of points to be placed by default on all connectors in 
the combustion chamber. This last parameter should be set equal to the number of points placed by 
default on connectors in the bottom of the swirler. 
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Procedure Calls 


This section calls the utility, large-scale creation, mesh creation, mesh refinement, and boundary 
condition procedures in such an order as to create a complete slice of the radial swirler. A slice is created 
to lessen the time needed to devise a solution for this case. One slice follows a single blade. The order of 
these procedure calls should not be modified. Some procedures rely upon other procedures already having 
been called. For example, the procedures CreateBuffer and CreateCombustionChamber assume that the 
procedure CreateSwirler has already been successfully called, and they will fail if it has not. 

Mesh Examples 

Figures 4 to 7 show various meshes that have been generated by the BladeMaker script. The 
parameters used to generate the meshes are typical for a LDI CFD design study. All the meshes comprise 
about 125 000 computational cells and, in general, take about 10 sec to generate. Using meshes with 
higher cell counts can take up to 2 min to generate. A Linux workstation with 2-GFlz Intel Xeon 5100 
series dual-core processors was used for mesh generation. 



Figure 4. — Top-down view of BladeMaker-generated grid with no blade tilt. 
Number of blades, 8; blade width, 0.2. 
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Figure 5. — Top-down view of BladeMaker-generated grid with blade tilt of 
1 5.0°. Number of blades, 8; blade width, 0.2. 



Figure 6. — Perspective view of BladeMaker-generated grid with blade tilt of 15.0°. Number 
of blades, 8; blade width, 0.2; swirler height, 1.0; combustion chamber height, 30.0. 
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Figure 7. — Side view of BladeMaker generated-grid with blade tilt of 1 5.0°. Number of blades, 8; blade width, 0.2; 
swirler height, 3.0; combustion chamber height, 30.0. 


Concluding Remarks 

The BladeMaker script removes a major time bottleneck for studying lean- direct-injection (LDI) 
concepts with computational fluid dynamics (CFD) codes, like the National Combustion Code (NCC). 
Although the geometry produced is a simple can-annular (canned) combustor, it is suitable for studying 
important combustor design parameters, such as the swirler blade angle. The canned combustor may be 
studied as a periodic slice or as a full can-annular combustor. The mesh density may also be varied, 
which is also important for determining whether or not the resulting CFD solution is mesh independent. 

We believe that this script is mostly error free in its current format, but 100-percent reliability cannot 
be guaranteed. Do not use this script in its current fonn for mission-critical applications. The script is 
being released in an open source licensing model, as documented in NASA OPEN SOURCE 
AGREEMENT VERSION 1.3.2 (ref. 16). By releasing BladeMaker in an open source model, we 
strongly believe that this will increase software quality via community peer review, accelerate 
BladeMaker development via community contributions, maximize the awareness and impact of NASA 
research, and increase dissemination of this software in support of NASA's education mission. 
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Appendix — BladeMaker Code Listing 


################################################################## 

################################################################## 


#### BLADEMAKER #### 

#### #### 

#### #### 

#### BladeMaker is a Glyph macro for Gridgen Ver 15.07 to be #### 
#### used for the automatic creation of radial swirlers. #### 

#### #### 

#### Scripting By: Daniel Thompson #### 

#### Grid & Geometry By: Anthony C. Iannetti #### 

#### Build Info: Release Version 1.0.2 #### 

#### Created: 18 July 2005 #### 

#### #### 


#### See manual for contact info. #### 

################################################################## 
################################################################## 

# Require the appropriate Glyph version. 

# Removal of this line may allow BladeMaker to run on Gridgen vl5.05 - 15.06. 
package require PWI Glyph 1.6.7 

######################################## 

# PARAMETERS # 

######################################## 

# Set Pi - Tel uses radians in trig functions, 
set PI [expr {4.0 * atan (1.0)}] 

# POSSIBLE VALUES FOR RUN_MODE 

# 

# 0: BladeMaker takes parameter values from params . in . Be sure 

# to use this mode if running this macro from unix-unopt . sh 

# 1: BladeMaker uses a crude GUI to get user input, 
set RUN_MODE 0 

set PARAM FILE "/your/working/directory/param. in" 
set EOF "^EOF*" 

# Read the parameter file, 
set cin [open $PARAM FILE] 

gets $cin line 

while {$line != $EOF}\ 

{ 

switch $line\ 

{ 

"ASW TYPE" {set $line [gets $cin] } 

"ASW DIM" {set $line [gets $cin] } 

"swirler rad" {set $line [gets $cin] } 

"bladevol_rad" {set $line [gets $cin] } 

"manif old_rad" {set $line [gets $cin] } 

"num_blades" {set $line [gets $cin] } 
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"blade_ang" {set $line [gets $cin] } 

"blade_rad" {set $line [gets $cin] } 

"height_s" {set $line [gets $cin] } 

"basic_con_dim" {set $line [gets $cin] } 

"blade_con_dim" {set $line [gets $cin] } 

"wall_ds" {set $line [gets $cin] } 

"height_cc" {set $line [gets $cin] } 

"fname" {set $line [gets $cin] } 

} 

gets $cin line 

} 

######################################## 

# UTILITY PROCEDURES # 

######################################## 

proc ClearWorkspace { } \ 

{ 

############################################## 

# Procedure: ClearWorkspace # 

# Purpose: Prepares the workspace for our # 

# use. # 

############################################## 

global ASW_TYPE 
global ASWJDIM 

gg : : memClear 

gg : : aswDeleteBC -glob "*" 

gg : : aswDeleteVC -glob "*" 

gg : : aswSet $ASW_TYPE -dim $ASW_DIM 

gg : : defReset 

gg : : tolReset 

gg: : updatePolicy DISPLAY_AND_INPUT 

} 

######################################## 

# CONNECTOR PROCEDURES # 

######################################## 

proc SplitConlntoNEqualParts {num pieces con list}\ 

{ 

############################################## 

# Procedure: SplitConlntoNEqualParts # 

# Purpose: Splits up connectors into equal # 

# parts. # 

# Note: This procedure adapted from # 

# ConSplitlntoN on the Pointwise Glyph # 

# exchange @ http://www.pointwise.com. # 

############################################## 

set return con list $con list 

foreach con $con_list\ 

{ 

for {set split 1} {$split < $num_pieces} finer split}\ 
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{ 

set ds [expr 1.0 / [expr {$num_pieces - $split + 1}]] 

set pt [gg : : conGetPt $con -arc $ds] 

set new_con [gg : : conSplit Scon $pt] 

set new dim [gg::conDim $new con] 

set needed dim [expr [$num pieces - $split + 1}] 

if {$new dim > 0 && $new dim < $needed dim}\ 

{ 

set delta [expr [$needed_dim - $new_dim}] 
set old dim [gg::conDim Scon] 

gg : : conRedimBegin 

gg::conRedim $new con [expr [$new dim + $delta}] 
gg::conRedim Scon [expr [$old dim - $delta}] 
gg : : conRedimEnd 
} 

set con $new con 

set return con list [lappend return^con list Scon] 

} 

} 

return Sreturn con list 

} 

proc CreateEndptCon { index }\ 

{ 

############################################## 


# Procedure: CreateEndptCon # 

# Purpose: Used to create a connector # 

# between the endpoints of two # 

# other connectors . # 


############################################## 

for [set i 0} [$i < [llength $index] } finer i 2}\ 

{ 

gg : : conBegin 
gg : : segBegin 

gg : : segAddControlPt [gg :: conGetPt [lindex [gg : : conGetAll ] [lindex $index $i] ] 
-arc 1.0] 

gg :: segAddControlPt [gg :: conGetPt [lindex [gg :: conGetAll ] [lindex $index 
[expr [$i + 1}]]] -arc 1.0] 
gg : : segEnd 

lappend r [gg::conEnd] 

} 

return $r 

} 

proc Create2PtCon [points }\ 

{ 

############################################## 

# Procedure: Create2PtCon # 

# Purpose: Used to create a connector # 

# between any 2 given points. # 

############################################## 
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for {set i 0} {$i < [llength $points]} (incr i 2}\ 

{ 

gg : : conBegin 
gg : : segBegin 

gg : : segAddControlPt [lindex $points $i] 

gg : : segAddControlPt [lindex $points [expr {$i + 1}]] 

gg : : segEnd 

lappend r [gg::conEnd] 

} 

return $r 

} 

proc SplitConAtPer {percent index} \ 

{ 

############################################## 

# Procedure: SplitConAtPer # 

# Purpose: Used to split a connector in half # 

# at a given percentage of its # 

# total length. # 

############################################## 

gg::conSplit [lindex [gg : : conGetAll ] $index] [gg : : conGetPt [lindex 

[gg : : conGetAll ] $index] -arc $percent] 

} 

######################################## 

# DATABASE PROCEDURES # 

######################################## 

proc Create2PtDB { index} \ 

{ 

############################################## 


# Procedure: Create2PtDB # 

# Purpose: Used to create a 2 point, # 

# straight, 3D line as a database # 

# entity. # 


############################################## 

for {set i 0} {$i < [llength $index] } {incr i 2}\ 

{ 

gg : : dbCurveBegin -type AKIMA 

gg : : dbCurveAddPt -db [list 0 0 [lindex [gg : : dbGetAll ] [lindex $index $i] ] ] 
gg : : dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] [lindex $index [expr 

{ $i + 1}]]]] 

lappend r [gg : : dbCurveEnd] 

} 

return $r 

} 

######################################## 

# MESHING PROCEDURES # 

######################################## 

proc CreateStrucDom { index} \ 
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{ 

############################################## 

# Procedure: CreateStrucDom # 

# Purpose: Used to create a structured # 

# domain . # 

############################################## 

for {set i 0} {$i < [llength $index] } finer i 4}\ 

{ 

gg::domBegin -type STRUCTURED 
gg : : edgeBegin 

gg : : edgeAddCon [lindex [gg : : conGetAll ] [lindex $index $ i ] ] 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon [lindex [gg :: conGetAll ] [lindex $index [expr {$i + 1}]]] 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon [lindex [gg :: conGetAll ] [lindex $index [expr {$i + 2}]]] 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon [lindex [gg :: conGetAll ] [lindex $index [expr {$i + 3}]]] 

gg : : edgeEnd 

lappend r [gg::domEnd] 

} 

return $r 

} 

proc CreateStrucBlk { index} \ 

{ 

############################################## 

# Procedure: CreateStrucBlk # 

# Purpose: Used to create a structured block.# 

############################################## 

for [set i 0} {$i < [llength $index] } finer i 6}\ 

f 


gg 

: blkBegin -type STRUCTURED 







gg 

: f aceBegin 










gg 

: f aceAddDom 

[lindex 

[gg: 

: domGetAll ] 

[lindex 

$index 

$i] ] 




gg 

: f aceEnd 










gg 

: f aceBegin 










gg 

: f aceAddDom 

[lindex 

[gg: 

: domGetAll ] 

[lindex 

$index 

[expr 

f $i 

+ 

1}] ] ] 

gg 

: f aceEnd 










gg 

: f aceBegin 










gg 

: f aceAddDom 

[lindex 

[gg: 

: domGetAll ] 

[lindex 

$index 

[expr 

f $i 

+ 

2}] ] ] 

gg 

: f aceEnd 










gg 

: f aceBegin 










gg 

: f aceAddDom 

[lindex 

[gg: 

: domGetAll ] 

[lindex 

$index 

[expr 

f$i 

+ 

3}] ] ] 

gg 

: f aceEnd 










gg 

: f aceBegin 










gg 

: f aceAddDom 

[lindex 

[gg: 

: domGetAll ] 

[lindex 

$index 

[expr 

f$i 

+ 

4}] ] ] 

gg 

: f aceEnd 










gg 

: f aceBegin 










gg 

: f aceAddDom 

[lindex 

[gg: 

: domGetAll ] 

[lindex 

$index 

[expr 

f $i 

+ 

5}] ] ] 

gg 

: f aceEnd 











lappend r [gg::blkEnd] 
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} 


return $r 

} 

######################################## 

# BOUNDARY CONDITION PROCEDURES # 

######################################## 

proc TagBoundaryConditions { } \ 

{ 

############################################## 

# Procedure: TagBoundaryConditions # 

# Purpose: Tags the boundary conditions. # 

############################################## 

### 

### Inlet-UVW (Bottom of Blades) 

### 


set doms [lrange [gg : : domGetAll ] 0 16] 
lappend doms [lindex [gg :: domGetAll ] 81] 
lappend doms [lindex [gg :: domGetAll ] 82] 
gg::aswSetBC $doms "inlet-uvw" 
unset doms 


### 

### Exit (Top of Combustion Chamber) 

### 


set doms [lrange [gg: 

: domGetAll] 144 

147] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

110] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

115] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

119] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

122] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

126] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

129] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

132] 


gg::aswSetBC $doms "exit" 
unset doms 


### 

### Periodic-A (Left Hand Swirler Cup Face) 

### 


lappend doms 
lappend doms 
lappend doms 
lappend doms 
lappend doms 
lappend doms 
gg : : aswSetBC 
unset doms 


[lindex [gg :: domGetAll ] 
[lindex [gg :: domGetAll ] 
[lindex [gg :: domGetAll ] 
[lindex [gg :: domGetAll ] 
[lindex [gg :: domGetAll ] 
[lindex [gg :: domGetAll ] 
$doms "periodic-a" 


34] 

78] 

91] 

105] 

116] 

133] 


### 

### Periodic-B (Right Hand Swirler Cup Face) 

### 
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lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

35] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

80] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

102] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

108] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

112] 

lappend 

doms 

[lindex 

[gg 

: domGetAll ] 

124] 


gg::aswSetBC $doms "periodic-b" 
unset doms 


### 

### Periodic-C (Blade Reigon Left Face) 

### 

gg : : aswCreateBC "BndCond 20" -solid 1 -id 20 
lappend doms [lindex [gg : : domGetAll ] 46] 
lappend doms [lindex [gg :: domGetAll ] 64] 
gg::aswSetBC $doms "BndCond 20" 
unset doms 

### 

### Periodic-D (Combustion Chamber Left Face) 

### 

gg :: aswCreateBC "BndCond 21" -solid 1 -id 21 
gg::aswSetBC [lindex [gg :: domGetAll ] 135] "BndCond 21 

### 

### Periodic-E (Blade Reigon Right Face) 

### 

gg :: aswCreateBC "BndCond 22" -solid 1 -id 22 
lappend doms [lindex [gg :: domGetAll ] 51] 
lappend doms [lindex [gg :: domGetAll ] 68] 
gg::aswSetBC $doms "BndCond 22" 

### 

### Periodic-F (Combustion Chamber Right Face) 

### 

gg :: aswCreateBC "BndCond 23" -solid 1 -id 23 
gg::aswSetBC [lindex [gg :: domGetAll ] 139] "BndCond 23 
} 

proc ExportASW [fname}\ 

{ 

############################################## 

# Procedure: ExportASW # 

# Purpose: Creates an analysis software # 

# file (like patran.out) . # 

############################################## 

gg : : aswExport $fname 

} 


######################################## 
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# MESH REFINEMENT PROCEDURES # 

######################################## 

proc ClusterPts {ds index loc}\ 

{ 

############################################## 

# Procedure: ClusterPts # 

# Purpose: Used to cluster points on the # 

# start or end of given connectors. # 

############################################## 

foreach i $index\ 

{ 

lappend pts [lindex [gg : : conGetAll ] $i] 

} 

if {$loc == "start" }\ 

{ 

gg : : conBeginSpacing $pts $ds 

} 

if {$loc == "end"}\ 

{ 

gg : : conEndSpacing $pts $ds 

} 

} 

proc Ref ineSwirlerGrid {wall ds blade con dim}\ 

{ 

############################################## 


# Procedure: Ref ineSwirlerGrid # 

# Purpose: Used to refine the grid of the # 

# swirler portion of the radial # 

# swirler. # 


############################################## 

# Redimension the reigon around the blades. 


gg: 

gg: 

: conRedimBegin 
: conRedim [lindex 

[gg: 

: conGetAll ] 

0] 

$blade con dim 

gg: 

:conRedim [lindex 

[gg: 

: conGetAll ] 

2] 

$blade con dim 

gg: 

: conRedimEnd 






# Cluster points around the blade wall. 

ClusterPts $wall_ds [list 31 33 34 64 68 73] "start" 
ClusterPts $wall_ds [list 30 35 36 62 65 70] "end" 

# Cluster points around the rear face, 
for [set i 88} [$i <= 115} finer i}\ 

{ 

lappend index $i 

} 

ClusterPts $wall_ds $index "start" 

# Cluster points around the forward face. 

ClusterPts $wall ds $index "end" 

unset index 
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} 


proc Ref ineCombustionChamberGrid {wall ds } \ 

{ 

############################################## 

# Procedure: Ref ineCombustionChamberGrid # 

# Purpose: Used to refine the grid of the # 

# combustion chamber portion of the # 

# radial swirler. # 

############################################## 

# Cluster points around the wall facing the blades, 
for {set i 156} {$i <= 159} finer i}\ 

{ 

lappend subcon index $i 

} 

for {set i 163} {$i <= 164} finer i}\ 

{ 

lappend subcon index $i 

} 

for {set i 171} {$i <= 172} finer i}\ 

{ 

lappend subcon index $i 

} 

for {set i 181} {$i <= 184} finer i}\ 

{ 

lappend subcon index $i 

} 

for {set i 193} {$i <= 194} finer i}\ 

{ 

lappend subcon index $i 

} 

lappend subcon index 175 
lappend subcon index 178 
lappend subcon index 167 

ClusterPts $wall ds $subcon index "start" 
unset subcon index 

# Cluster points around the wall near the outside of the chamber, 
for {set i 195} {$i <= 204} {incr i}\ 

{ 

lappend subcon index $i 

} 

ClusterPts $wall ds $subcon index "end" 

unset subcon index 

} 

# Ref ineBuf f erGrid taken out - 13 July 2005 
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# CreateBuf f erGrid now assumes that Ref ineSwirlerGrid 

# has already been called, thus CreateBuf ferGrid 

# will copy the refinement. 

######################################## 

# LARGE-SCALE CREATION PROCEDURES # 

######################################## 

proc MakeWhole {num blades }\ 

{ 

############################################## 

# Procedure: MakeWhole # 

# Purpose: Copies the geometry and grid via # 

# rotation about the z-axis to # 

# generate the whole radial swirler.# 

############################################## 

set db xform [gg : : dbGetAll ] 

set blk xform [gg : : blkGetAll ] 

set wedge_ang [expr {360.0 / $num_blades } ] 

for {set i 1} {$i < $num_blades} finer i}\ 

{ 

gg : : dbCopyBegin $db xform 

gg : : xf ormRotate [list 000] [list 001] [expr [$wedge_ang * $i}] 
gg : : dbCopyEnd 

gg : : blkCopyBegin $blk xform 

gg :: xf ormRotate [list 000] [list 001] [expr { $wedge_ang * $i}] 
gg : :blkCopyEnd 
} 

} 

proc CreateSwirler [swirler_rad bladevol rad manifold rad num blades 
blade ang blade rad height basic con dim}\ 

{ 

############################################## 

# Procedure: CreateSwirler # 

# Purpose: Creates the swirler cup-blade # 

# region of the radial swirler. # 

############################################## 

# VARIABLES 

# Miscellaneous 
global PI 

global VERBOSEJMODE 
global RUN_MODE 

set wedge_ang [expr {(2.0 * $PI) / $num_blades } ] 

# Lower Plane Points 


set origin [list 

0 0 

0] 








set swirler lp [ 

[list 

[expr {-$swirler 

rad 

* 

(sin 

($wedge ang 

/ 2.0) 

) }] 

[expr 

[$swirler rad * 

(cos 

($wedge ang / 2.0 

) ) }] 

0] 






set swirler cp [ 

list 

0 $swirler rad 0] 








set swirler rp 

[list 

[expr [$swirler 

rad 

★ 

(sin 

($wedge ang 

/ 2.0) 

) }] 

[expr 

[$swirler rad * 

(cos 

($wedge ang / 2.0 

) ) }] 

0] 
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set bladevol_lp [list [expr { -$bladevol_rad * (sin ($wedge_ang / 2.0))}] 

[expr { $bladevol_rad * (cos ($wedge_ang / 2.0))}] 0] 

set bladevol_cp [list 0 $bladevol_rad 0] 

set bladevol_rp [list [expr { $bladevol_rad * (sin ($wedge_ang / 2.0))}] [expr 

{ $bladevol_rad * (cos ($wedge_ang / 2.0))}] 0] 

set manifold lp [list [expr {-$manifold rad * (sin ($wedge ang / 2.0))}] 

[expr { $manif old_rad * (cos ($wedge_ang / 2.0))}] 0] 

set manifold cp [list 0 $manifold rad 0] 

set manifold rp [list [expr [$manifold rad * (sin ($wedge ang / 2.0))}] [expr 

{ $manif old_rad * (cos ($wedge_ang / 2.0))}] 0] 

set inner blade lp [list -$blade rad [expr [$swirler rad + 2 * $blade rad}] 

0] 

set inner_blade_cp [list 0 [expr { $swirler_rad + $blade_rad}] 0] 

set inner blade rp [list $blade rad [lindex $inner blade_lp 1] 0] 

set outer blade lp [list [lindex $inner blade lp 0] [expr [$bladevol rad - 2 

# $blade_rad}] 0] 

set outer_blade_cp [list 0 [expr { $bladevol_rad - $blade_rad}] 0] 

set outer blade rp [list [lindex $inner blade rp 0] [lindex $outer blade lp 

1] 0] 

# LOWER PLANE 

# Define Database Entities 


# Add Points to 
gg : rdbPtsBegin 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : rdbPtsAddPt 
gg : : dbPtsEnd 


Database 

$origin 
$swirler lp 
$swirler cp 
$swirler rp 
$bladevol lp 
$bladevol_cp 
$bladevol rp 
$manifold lp 
$manifold cp 
$manifold rp 
$inner blade lp 
$inner blade cp 
$inner blade rp 
$outer_blade_lp 
$outer_blade_cp 
$outer_blade_rp 


Create2PtDB [list 0 1\ 

1 4\ 

4 7 \ 

0 3\ 

3 6\ 

6 9] 


# Swirler Domain 

- Arc 








gg 

: dbCurveBegin 

-type 

o 
1— 1 

o 

a 

t- 1 

> 

ARC 




gg 

: dbCurveAddPt 

-db [ 

list 

0 

0 

] lindex 

[gg: 

: dbGetAll ] 

1] ] 

gg 

: dbCurveAddPt 

-db [ 

list 

0 

0 

[ lindex 

[gg: 

: dbGetAll ] 

3] ] 

gg 

gg 

: dbCurveAddPt 
: dbCurveEnd 

-db [ 

list 

0 

0 

] lindex 

[gg: 

: dbGetAll ] 

2] ] 
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Create2PtDB [list 10 13\ 

12 15] 

# Blade - Innermost Arc 

gg : : dbCurveBegin -type CIRCULAR ARC 

gg : : dbCurveAddPt -db [list 0 0 [lindex [gg : : dbGetAll ] 10]] 

gg : : dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 12]] 

gg :: dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 11]] 

gg : rdbCurveEnd 

# Blade - Outermost Arc 

gg :: dbCurveBegin -type CIRCULAR ARC 

gg :: dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 13]] 

gg :: dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 15]] 

gg :: dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 14]] 

gg : rdbCurveEnd 

# Rotate Left Blade Volume and Left Manifold Volume Akimas 

gg : : dbTransf ormBegin [list [lindex [gg :: dbGetAll ] 17] [lindex [gg :: dbGetAll ] 

18] [lindex [gg :: dbGetAll ] 4] [lindex [gg :: dbGetAll ] 7]] -maintain linkage 

gg : : xf ormRotate $swirler lp [ggu : : vec3Add $swirler lp [list 0 0 1]] 

$blade_ang 

gg : : dbTransf ormEnd 

# Rotate Right Blade Volume and Right Manifold Volume Akimas 

gg :: dbTransf ormBegin [list [lindex [gg :: dbGetAll ] 20] [lindex [gg :: dbGetAll ] 

21] [lindex [gg :: dbGetAll ] 6] [lindex [gg :: dbGetAll ] 9]] -maintain linkage 

gg :: xf ormRotate $swirler rp [ggu : : vec3Add $swirler_rp [list 0 0 1]] 

$blade_ang 

gg : : dbTransf ormEnd 

# Rotate Blade Volume Center Point 

gg :: dbTransf ormBegin [list [lindex [gg :: dbGetAll ] 5] [lindex [gg :: dbGetAll ] 

8] ] -maintain linkage 

gg :: xf ormRotate $swirler cp [ggu : : vec3Add $swirler cp [list 0 0 1]] 

$blade_ang 

gg : : dbTransf ormEnd 

# Rotate Blade 

gg :: dbTransf ormBegin [list [lindex [gg :: dbGetAll ] 23] [lindex [gg :: dbGetAll ] 

24] [lindex [gg :: dbGetAll ] 26] [lindex [gg :: dbGetAll ] 25] [lindex 

[gg :: dbGetAll ] 13] [lindex [gg :: dbGetAll ] 14] [lindex [gg :: dbGetAll ] 15] 

[lindex [gg :: dbGetAll ] 10] [lindex [gg :: dbGetAll ] 12]] -maintain linkage 

gg :: xf ormRotate $inner blade cp [ggu : : vec3Add $inner blade cp [list 001]] 

$blade_ang 

gg : : dbTransf ormEnd 

# Blade Volume - Arc 

gg :: dbCurveBegin -type CIRCULAR ARC 

gg :: dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 4]] 

gg :: dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 6]] 

gg :: dbCurveAddPt -db [list 0 0 [lindex [gg :: dbGetAll ] 5]] 

gg : rdbCurveEnd 

# Manifold Volume - Arc 

gg :: dbCurveBegin -type CIRCULAR ARC 
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gg: 

: dbCurveAddPt 

-db 

gg: 

: dbCurveAddPt 

-db 

gg: 

: dbCurveAddPt 

-db 

gg: 

: dbCurveEnd 



[list 0 0 [lindex 
[list 0 0 [lindex 
[list 0 0 [lindex 


[gg : : dbGetAll ] 7]] 
[gg : : dbGetAll ] 9]] 
[gg :: dbGetAll ] 8]] 


# Define Connectors 

# Create the connectors on the database entities, 
gg : : conOnDBEnt ALL 


# Swirler Cup 

SplitConlntoNEqualParts 4 [list [lindex [gg : : conGetAll ] 6]] 
SplitConAtPer 0.25 0 
SplitConAtPer 0.25 2 

Create2PtCon [list $origin $swirler_cp] 

SplitConAtPer 0.4 18 
CreateEndptCon [list 14 18\ 

16 18] 

gg : : conDelete [lindex [gg :: conGetAll ] 18] 

SplitConAtPer 0.5 19 
SplitConAtPer 0.5 19 
SplitConAtPer 0.5 14 
SplitConAtPer 0.5 15 
CreateEndptCon [list 21 19\ 

23 17] 

SplitConAtPer 0.5 25 
SplitConAtPer 0.5 25 
CreateEndptCon [list 17 10\ 

19 12] 


4 [list [lindex [gg :: conGetAll ] 8]] 
2 [lrange [gg :: conGetAll ] 6 7] 

28\ 


# Blade Volume and Blade 
SplitConlntoNEqualParts 
SplitConlntoNEqualParts 
CreateEndptCon [list 7 

8 32\ 

34 2 9\ 

9 30\ 

4 36\ 

5 39\ 

9 33] 

Create2PtCon [list [gg : : conGetPt [lindex [gg :: conGetAll ] 

[gg : : conGetPt [lindex [gg :: conGetAll ] 32] -arc 0.0]] 


# Manifold Volume 

SplitConlntoNEqualParts 4 [lindex [gg :: conGetAll ] 6] 

for [set i 27} [$i <= 29} finer i}\ 

{ 

CreateEndptCon [list $i [expr [16 + $i}]] 

} 

gg::conDim ALL $basic con dim 

# Create Domains 

# Swirler Cup 
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CreateStrucDom [list 19 23 21 17\ 


23 

20 

15 

22\ 

22 

16 

14 

24\ 

21 

24 

13 

18\ 

13 

25 

6 

10\ 

14 

12 

7 

25\ 

16 

26 

8 

12\ 

15 

11 

9 

26] 


# Blade Volume and Blade 
CreateStrucDom [list 6 35 27 0\ 


42 

4 

39 

35\ 

7 

36 

31 

42\ 

39 

33 

37 

28\ 

8 

41 

32 

36\ 

34 

40 

29 

37\ 

41 

38 

40 

5\ 

9 

2 

30 

38] 


# Manifold 

CreateStrucDom [list 27 47 43 1\ 


28 

48 

44 

47\ 

29 

49 

45 

4 8\ 

30 

3 

46 

49] 


# UPPER PLANE 

# Create Database Entities via Copy 
gg : : dbCopyBegin ALL 

gg : : xf ormTranslate "0 0 $height" 
gg : : dbCopyEnd 

# Create Connectors and Domains via Copy 
gg : : domCopyBegin ALL 

gg :: xf ormTranslate "0 0 $height" 
gg : rdomCopyEnd 

# MIDDLE PLANE 

# Create Database Entities and Connectors 

gg::conDim [gg : : conOnDBEnt [Create2PtDB [list 0 29\ 

1 30\ 

2 31\ 

3 32\ 

4 33\ 

5 34 \ 

6 35\ 

7 36\ 

8 37\ 

9 38] ] ] $basic_con_dim 

gg::conDim [CreateEndptCon [list 17 53\ 

19 50\ 

21 52 \ 

18 61 \ 

20 54 \ 

13 60\ 
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22 56\ 

14 58\ 

26 67\ 

25 62 \ 

35 7 1 \ 

38 87\ 

42 74 \ 

32 83\ 

31 78\ 

4 7 5\ 

5 88\ 

33 7 9\ 

43 92\ 

45 97]] $basic con dim 


# Create Domains 
CreateStrucDom [list 17 100 53 


18 

110 

61 

113\ 

19 

111 

50 

100\ 

10 

113 

64 

101\ 

114 

69 

103 

11\ 

6 

119 

63 

101\ 

7 

102 

66 

1 1 9\ 

8 

118 

68 

102\ 

9 

103 

70 

1 1 8\ 

23 

111 

51 

112\ 

24 

112 

59 

115\ 

15 

114 

55 

1 1 6\ 

116 

57 

117 

1 6\ 

25 

115 

62 

1 1 9\ 

12 

117 

65 

102\ 

26 

116 

67 

118\ 

0 

101 

73 

1 04 \ 

119 

71 

120 

35\ 

4 

122 

75 

125\ 

5 

123 

88 

12 6\ 

38 

118 

87 

12 1 \ 

103 

89 

106 

2\ 

42 

119 

74 

122\ 

31 

124 

78 

122\ 

124 

83 

123 

32\ 

118 

82 

123 

41\ 

27 

104 

72 

120\ 

28 

120 

81 

1 05\ 

105 

86 

121 

2 9\ 

120 

76 

125 

39\ 

125 

79 

127 

33\ 

127 

34 

126 

84\ 

85 

121 

40 

12 6\ 

121 

90 

106 

30\ 

104 

93 

107 

1\ 

120 

91 

128 

47\ 

105 

94 

108 

48\ 

121 

96 

129 

4 9\ 

106 

98 

109 

3\ 

107 

92 

128 

43\ 

128 

95 

108 

44\ 


110 \ 
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108 97 129 45\ 

129 99 109 4 6\ 

127 80 105 37\ 

102 77 124 36\ 

21 110 52 112\ 

112 56 116 22 \ 

20 111 54 114\ 

13 113 60 115\ 

115 58 117 14] 

# CREATE ALL STRUCTURED BLOCKS 

CreateStrucBlk [list 0 20 40 42 49 85\ 


3 

23 

41 

50 

85 

88\ 

1 

21 

49 

51 

86 

87\ 

2 

22 

50 

52 

86 

89\ 

4 

24 

43 

45 

53 

8 8 \ 

5 

25 

46 

53 

54 

8 9\ 

6 

26 

47 

52 

54 

55\ 

7 

27 

44 

48 

51 

55\ 

8 

28 

45 

56 

57 

66\ 

9 

29 

57 

58 

62 

69\ 

10 

30 

46 

62 

63 

84 \ 

11 

31 

67 

69 

70 

83\ 

13 

33 

68 

71 

72 

83\ 

12 

32 

47 

64 

65 

84 \ 

14 

34 

59 

60 

65 

72\ 

15 

35 

48 

60 

61 

7 3\ 

16 

36 

66 

74 

75 

7 9\ 

17 

37 

67 

75 

76 

80\ 

18 

38 

68 

76 

77 

8 1 \ 

19 

39 

73 

77 

78 

82] 


### 

### BEGIN PATCH 

### 

### 

### 

gg : : blkDelete [lrange [gg : : blkGetAll ] 0 2] -doms -cons 
gg : : conOnDBEnt [lindex [gg : : dbGetAll ] 58] 

gg : : conBegin 
gg : : segBegin 

gg : : segAddControlPt $origin 

gg : : segAddControlPt [gg : : conGetPt [lindex [gg : : conGetAll ] 98] -arc 0.0] 

gg : : segEnd 
gg : : conEnd 

gg : : conBegin 
gg : : segBegin 

gg :: segAddControlPt $origin 

gg :: segAddControlPt [gg :: conGetPt [lindex [gg :: conGetAll ] 99] -arc 0.0] 
gg : : segEnd 
gg : : conEnd 

gg : : conBegin 
gg : : segBegin 


- ADDED 13 JULY 2005 TO AVOID HAVING A CELL 
AGAINST 2 PERIOD BOUNDARY CONDITIONS 
PATCH IS HACKY - REMOVE IN FUTURE 
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gg : : segAddControlPt $origin 

gg : : segAddControlPt [gg : : conGetPt [lindex [gg : : conGetAll ] 18] -arc 0.0] 

gg : : segEnd 
gg : : conEnd 

CreateEndptCon [list 115 98 \ 

115 97 \ 

115 99 ] 

gg::conDim [lrange [gg :: conGetAll ] 115 121] $basic con dim 

CreateStrucDom [list 115 119 98 116\ 

115 120 97 118\ 

115 121 99 117\ 

116 118 18 13 \ 

118 117 15 17 \ 

119 49 48 120\ 

120 121 44 45 ] 

CreateStrucBlk [list 40 76 78 79 81 83\ 

41 75 79 80 82 84] 

### 

### END SWIRLER GRID PATCH 

### 

} 

proc CreateBuffer [height of swirler cup}\ 

{ 

############################################## 


# Procedure: CreateBuffer # 

# Purpose: Creates the buffer between the # 

# combustion chamber and the # 

# swirler cup. # 


############################################## 

gg : : dbCopyBegin [concat\ 

[lrange [gg : : dbGetAll ] 29 32] \ 

[lindex [gg :: dbGetAll ] 45] \ 

[lindex [gg :: dbGetAll ] 48] \ 

[lindex [gg :: dbGetAll ] 51] \ 

[lrange [gg :: dbGetAll ] 58 61]] 

gg : : xf ormTranslate "0 0 $height_of swirler cup" 
gg : : dbCopyEnd 

gg : : blkCopyBegin [concat [lrange [gg : : blkGetAll ] 0 4] [lrange [gg : : blkGetAll ] 

17 18]] 

gg :: xf ormTranslate "0 0 $height_of swirler cup" 
gg : :blkCopyEnd 
} 

proc CreateCombustionChamber [height of_swirler cup 

height of combustion chamber basic con dim}\ 

{ 

############################################## 

# Procedure: CreateCombustionChamber # 
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# Purpose: Creates the combustion chamber. # 

# Note: Assumes CreateSwirler has been # 

# called. # 

############################################## 

set twice height of swirler cup [expr [2 * $height of swirler cup}] 

# Figure out the factor to multiply basic con dim 

# by to maintain grid density. 

if {$height of swirler cup < $height of combustion chamber }\ 

{ 

set chamber_con dim [expr { int ( ( $height of combustion chamber / 

$height of swirler cup) * $basic con dim * 0.15) }] 

}\ 

\ 

elseif {$height of swirler cup == $height of combustion chamber }\ 

{ 

set chamber con dim $basic con dim 

}\ 

\ 

else\ 

{ 

set chamber con dim [expr [int($height of swirler cup / 

$height of combustion chamber} * $basic_con dim * 0.15)] 

} 

# INNER COMBUSTION CHAMBER 

# Copy what we can from the buffer. 

gg : : dbCopyBegin [lrange [gg : : dbGetAll ] 68 78] 

gg : : xf ormTranslate "0 0 $height of swirler cup" 
gg : : dbCopyEnd 

gg : : domCopyBegin [concat [lindex [gg : : domGetAll ] 85] [lindex [gg : : domGetAll ] 

90] [lindex [gg :: domGetAll ] 95] \ 

[lindex [gg :: domGetAll ] 98] [lindex [gg :: domGetAll ] 101] [lindex 

[gg :: domGetAll ] 107] \ 

[lindex [gg :: domGetAll ] 109]] 
gg :: xf ormTranslate "0 0 $height of swirler cup" 
gg : :domCopyEnd 

# Create connectors on database entities. 

gg : : conOnDBEnt [lrange [gg :: dbGetAll ] 86 89] 

# Create extra connectors. 

CreateEndptCon [list 132 158\ 

141 1 63\ 

149 1 67\ 

150 168\ 

151 1 69\ 

130 156\ 

122 152\ 

124 154] 

gg::conDim [lrange [gg : : conGetAll ] 170 181] $basic con dim 

# Create Domains 
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CreateStrucDom [list 170 167 176 149\ 


170 

168 

177 

150\ 


170 

169 

178 

151\ 


176 

156 

179 

130\ 


177 

155 

179 

125\ 


177 

152 

180 

122\ 


178 

164 

180 

143\ 


179 

154 

181 

124\ 


180 

153 

181 

123\ 


176 

159 

171 

133\ 


179 

157 

174 

131\ 


181 

160 

172 

137\ 


180 

162 

175 

140\ 


178 

165 

173 

144\ 


158 

174 

132 

171\ 


161 

172 

138 

174\ 


172 

163 

175 

141\ 


175 

166 

173 

145] 


# Create 

Blocks 


CreateStrucBlk [list 

109 116 

107 

115 

117 

118 120 

12 1 \ 

85 

110 

121 

122 124 

125\ 

101 

114 

123 

129 130 

134\ 

98 

113 

125 

128 129 

133\ 

95 

112 

124 

127 128 

132\ 

90 

111 

120 

126 127 

131] 


# Stretch the Database to Size 

gg : : dbTransf ormBegin [lrange [gg : : dbGetAll ] 79 89] 

gg : : xf ormStretch [list 0 0 [expr [2 * $height of swirler cup}]]\ 

[list 0 0 [expr [3 * $height_of_swirler_cup } ] ] \ 

[list 0 0 [expr { (3 * $height_of_swirler_cup) 

$height_of combustion chamber}]] 
gg : : dbTransf ormEnd 

# Stretch the Blocks to Size 

gg :: blkTransf ormBegin [lrange [gg : : blkGetAll ] 26 32] 

gg :: xf ormStretch [list 0 0 [expr [2 * $height of swirler cup}]]\ 

[list 0 0 [expr [3 * $height_of_swirler_cup } ] ] \ 

[list 0 0 [expr { (3 * $height_of_swirler_cup) 

$height of combustion chamber}]] 
gg : : blkTransf ormEnd 

# OUTER CHAMBER 

# Copy some DB entities from the blade reigon and create connectors on 

# them. 

gg : : dbCopyBegin [concat\ 

[lrange [gg :: dbGetAll ] 7 9]\ 

[lindex [gg :: dbGetAll ] 28] \ 

[lrange [gg :: dbGetAll ] 36 38] \ 

[lindex [gg :: dbGetAll ] 57] \ 

[lrange [gg :: dbGetAll ] 65 67]] 

gg : : xf ormTranslate "0 0 $twice height of swirler cup" 
gg : : conOnDBEnt [gg : : dbCopyEnd] 
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SplitConlntoNEqualParts 2 [lrange [gg : : conGetAll ] 185 188] 

gg::conDim [lrange [gg :: conGetAll ] 182 192] $basic con dim 

gg::conDim [CreateEndptCon [list 185 189\ 

187 191]] $basic_con_dim 

# Stretch to fit. 

gg : : dbTransf ormBegin [lrange [gg : : dbGetAll ] 90 100] 

gg : : xf ormStretch [list 0 0 [expr [2 * $height of swirler cup}]]\ 

[list 0 0 [expr [3 * $height_of_swirler_cup } ] ] \ 

[list 0 0 [expr { (3 * $height_of_swirler_cup) 

$height_of combustion chamber}]] 
gg : : dbTransf ormEnd 

gg :: conTransf ormBegin [lrange [gg :: conGetAll ] 182 194] 

gg :: xf ormStretch [list 0 0 [expr [2 * $height of swirler cup}]]\ 

[list 0 0 [expr [3 * $height_of_swirler_cup } ] ] \ 

[list 0 0 [expr { (3 * $height_of_swirler_cup) 

$height of combustion chamber}]] 
gg : : conTransf ormEnd 

# Create some more DB entities, and create dimensioned connectors on them. 
gg::conDim [gg : : conOnDBEnt [Create2PtDB [list 80 94\ 

81 95\ 

82 96\ 

69 90\ 

70 91 \ 

71 92]]] $basic_con_dim 

# Create extra dimensioned connectors. 
gg::conDim [CreateEndptCon [list 179 189\ 

173 191\ 

132 1 85\ 

141 187]] $basic con dim 

# Create Domains 

CreateStrucDom [list 181 195 182 198\ 


178 

201 

193 

203\ 

175 

196 

183 

1 99\ 

171 

202 

194 

204\ 

172 

197 

184 

200\ 

182 

189 

193 

1 85\ 

193 

190 

183 

1 8 6\ 

183 

191 

194 

1 87\ 

194 

192 

184 

188\ 

179 

195 

189 

201\ 

176 

201 

190 

1 96\ 

173 

196 

191 

202\ 

169 

202 

192 

1 97\ 

132 

203 

185 

1 98\ 

138 

199 

186 

203\ 

141 

204 

187 

1 99\ 

145 

200 

188 

204] 

# Create 

Blocks 
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CreateStrucBlk [list 134 135 136 140 144 148\ 

131 136 137 141 145 149\ 

128 137 138 142 146 150\ 

125 138 139 143 147 151] 

# Redimension grid to maintain grid density, 
gg : : conRedimBegin 

gg::conRedim [lindex [gg : : conGetAll ] 195] $chamber con dim 
gg::conRedim [lindex [gg :: conGetAll ] 10] $chamber con dim 

gg::conRedim [lindex [gg :: conGetAll ] 6] $chamber con dim 

gg::conRedim [lindex [gg :: conGetAll ] 7] $chamber con dim 

gg::conRedim [lindex [gg :: conGetAll ] 8] $chamber con dim 

gg::conRedim [lindex [gg :: conGetAll ] 156] $chamber con dim 

gg : :conRedimEnd 

} 

######################################## 

# GRAPHICAL INTERFACE PROCEDURES # 

######################################## 

proc DrawButtonFireEvent { databaseAttriblnputSwirlerCupRadius 

databaseAttriblnputBladeVolumeRadius databaseAttriblnputManif oldRadius 

dat abaseAtt rib I nputNumberOf Blades databaseAttriblnputBladeAngle 

databaseAttriblnputBladeRadius databaseAttriblnputHeightS 

databaseAttriblnputHeightCC meshAttriblnputDef aultConnectorDim} \ 

{ 

############################################## 

# Procedure: DrawButtonFireEvent # 

# Purpose: Handle the pressing of the DRAW # 

# button in the GUI . # 

############################################## 

ClearWorkspace 

CreateSwirler [ $databaseAttribInputSwirlerCupRadius get] 

[$ dat abaseAtt rib I nputBladeVolumeRadius get] 

[ $dat abaseAtt rib I nputManif oldRadius get] [ $ dat abaseAtt rib I nputNumberOf Blades 
get] [ $databaseAttribInputBladeAngle get] [ $databaseAttribInputBladeRadius 
get] [ $databaseAttribInputHeightS get] [ $meshAttribInputDef aultConnectorDim 
get] 

CreateBuffer [ $databaseAttribInputHeightS get] 

CreateCombustionChamber [ $databaseAttribInputHeightS get] 

[$databaseAttribInputHeightCC get] [ $meshAttribInputDef aultConnectorDim get] 
MakeWhole [ $databaseAttribInputNumberOfBlades get] 

} 

proc CreateGUI [}\ 

{ 

############################################## 


# Procedure: CreateGUI # 

# Purpose: Creates a very crude GUI for # 

# BladeMaker. # 

# Note: This procedure is for TESTING # 

# PURPOSES ONLY. The GUI is cryptic # 

# and probably won't even exist in # 

# BladeMaker 1.0 # 


############################################## 
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gg : : tkLoad 


wm title . "BladeMaker Testing GUI" 
set row 0 

foreach i { SwirlerCupRadius BladeVolumeRadius Manif oldRadius NumberOfBlades 
BladeAngle BladeRadius Heights HeightCC}\ 

{ 

label . databaseAttribLabel$i -text $i 

entry . databaseAttribInput$i 

lappend entrylist . databaseAttribInput$i 

grid . databaseAttribLabel$i -row $row -column 0 -sticky w 
grid . databaseAttribInput$i -row $row -column 1 

incr row 

} 

set row 0 

foreach i { Def aultConnectorDim} \ 

{ 

label ,meshAttribLabel$i -text $i 

entry . meshAttribInput$i 

lappend entrylist ,meshAttribInput$i 

grid ,meshAttribLabel$i -row $row -column 2 -sticky w 
grid . meshAttribInput$i -row $row -column 3 

incr row 

} 

set row 8 
set col 0 

frame .buttons 

button . buttons . draw -text "Draw" -command "DrawButtonFireEvent $entrylist" 
button . buttons . close -text "Close" -command "exit" 

pack . buttons . draw -side left 
pack . buttons . close -side right 

grid .buttons -row $row -column $col 

gg : : tkLoop 

} 

######################################## 

# BEGIN MAIN GLYPH SCRIPT # 

######################################## 

### 

### DO ***NOT*** modify this section unless 

### YOU KNOW WHAT YOU ARE DOING 

### ORDER OF FUNCTION CALLS ***DOES*** MATTER 
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### 

ClearWorkspace 

switch $RUN_MODE\ 

{ 

0 \ 

{ 

CreateSwirler $swirler rad $bladevol rad $manifold rad $num blades $blade ang 

$blade rad $height s $basic con dim 

Ref ineSwirlerGrid $wall ds $blade con dim 

CreateBuffer $height s 

CreateCombustionChamber $height s $height cc $basic con dim 
Ref ineCombustionChamberGrid $wall ds 
TagBoundaryConditions 
ExportASW $fname 
} 

\ 

1 \ 

{ 

CreateGUI 

} 

} 

### 

### END DO NOT MODIFY SECTION 

### 
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